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Abstract 

We calculate the leading order diphoton generalized distribution amplitudes by calculating 
the amplitude of the process 7*7 — > 77 in the low energy and high photon virtuality region at 
the Born order and in the leading logarithmic approximation. As in the case of the anomalous 
photon structure functions, the 77 generalized distribution amplitudes exhibit a characteristic 
In Q 2 behaviour and obey inhomogeneous QCD evolution equations. 



1 Introduction 



The photon is a much interesting object for QCD studies. Its pointlike coupling to quarks enables to 
calculate perturbatively part of its wave function. In the case of a virtual photon, this perturbative 
part is leading at large virtuality; a twist expansion generates non-leading components of the photon 
distribution amplitude pQ, from which the lowest order one is chiral-odd and proportional to the 
magnetic susceptibility of the vacuum. The study of the two photon state is kinematically richer and 
is thus a most welcome theoretical laboratory for the study of exclusive hard reactions. 

The parton content of the photon has been the subject of many studies since the seminal paper 
by Witten [2]. In a recent paper [3] we extended the notion of anomalous parton distribution in a 
photon to the case of generalized parton distributions (GPDs) used for the factorized description of 
the non-diagonal kinematics of deeply virtual Compton scattering (DVCS), 7* (9)7 — > 77, namely at 
large energy and small hadronic momentum transfer but large photon virtuality (Q 2 = —q 2 ). 

The two meson GDAs generalize [H [5] the concept of usual distribution amplitude of mesons and 
baryons introduced long time ago in the studies of exclusive processes [B] . They describe the coupling 
of a quark-antiquark (or gluon-gluon) pair to a pair of mesons, and are related by crossing to the 
meson GPDs. In the same way, the two photon generalized distribution amplitudes (GDAs) which 
describes the coupling of a quark antiquark (or gluon-gluon) pair to a pair of photons, are related by 
crossing to the photon GPDs. 

We thus study here the scattering amplitude of the 7* (9)7 — > 77 process in the near threshold 
kinematics, namely at small s and large —t ~ Q 2 , at large Q 2 in the leading order of the electro- 
magnetic coupling. This enables us to define and calculate perturbatively the Born approximation 
of the diphoton GDAs. This Born order contribution enters the QCD evolution equations as an 
inhomogeneous term and we exhibit the properties of the evolved solution. 

On the phcnomcnlogical side, we do not expect any study of this process to be feasible since the 
production of hadronic states that decay into a pair of photons - mostly pseudoscalar neutral mesons 
- will strongly dominate the counting rates of this purely electromagnetic process. 

The analysis presented in this paper is reminiscent of the perturbative calculation of the two p 
GDAs in terms of the single p DAs [7J, for which both incoming photon were chosen to be hard in 
order to justify this factorization of the GDAs. 

2 The 7*(<?)7(pi) — » l{q')l{p2) process in the threshold region 

Two photon production in Compton scattering on a photon target 

i*(q)i(q ) -> 7(pibb2) (1) 

involves, at leading order in a em , and zeroth order in as the six Feynman diagrams of Fig [I] with 
quarks in the loop. Note that we do not consider in this paper the pure QED content of this amplitude, 
involving a lepton in the loop. It can straightforwardly be obtained from the results below. 

For simplicity, we restrict in this paper to the threshold kinematics where W 2 = (pi + P2) 2 = 0. 
This simplifies greatly the tensorial structure of the amplitude while still preserving the richness of 
the skewedness (£) dependence of GDAs, but forbids any impact parameter interpretation of the 
GDAs [SJ. Our conventions for the kinematics are the following: 




Pi = (P , P2 = (p , C = 1 - C) 

where p and n are two light-cone Sudakov vectors and 2p ■ n — s. The momentum I in the quark loop 
is parametrized as 

J" = zp» + fin" + l T , (2) 

with I j, — — l 2 . The process involves one virtual and three real photons and its amplitude can be 
written as 

A = e^e&^T^, (3) 
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Figure 1: The Born order diagrams for 7* 7 — > 7 7 

where in our (forward) kinematics the four photon polarization vectors e(g), e'(q'), ei(pi) and £2(^2) 
are transverse with respect to the Sudakov vectors p and n. 
The tensorial decomposition of T^ al3 reads [8] 

T^(W = 0) = + I + 9^9^ ~ 9"t9t P ) + \ (g^itf - W 3 , (4) 

and it involves three scalar functions Wj, i — 1, 2, 3. 

The integration over / is performed as usual within the Sudakov representation, using 

dH = ^dzd(3 d 2 l T dl 2 . 



Before starting the explicit computation, we note that in order to interpret our result in terms of 
factorized quantities, we will keep our expressions unintegrated with respect to z. In all the following 
calculations, the mass of the quark will play the role of an infrared regulator. Let us briefly present 
our calculations for diagram A, and then give final results for the other diagrams. We get for 



dzdfidl 2 TrA x 



[{I — q) 2 — m 2 + ir)][(l — pi) 2 — m 2 + irj\(l 2 — m 2 + 



- q') 2 — m 2 + irj\ ' 
(5) 

where TrA\ — Tr[jj,{l — q + m)^(l — p\ — P2 + m)j^(l — p\ + 771)7^ (7 + m)] and Tr means the trace 
over spinorial indices. 
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Wc first integrate over (3 using the Cauchy theorem. The propagators induce poles in the complex 
/3-plane with values 



l 2 + m 2 — ir/ l 2 + m 2 — if] 

Pi = > P2 = — 7 tt , 

zs [z — Q)s 

l 2 + m 2 + zQ 2 - 177 n l 2 + m 2 - irj 

03 = = , P4 = Z ■ (6) 

zs zs 

Since the four poles lie all below the real axis for z > 1 and lie all above the real axis for z < 0, the 
only region where the amplitude may not vanish is 1 > z > 0. One identifies two different regions: 

• the region where £ < z < 1, for which one may close the contour in the lower half plane and get 
the contribution of the poles p\ and 02 , 

• the region where < z < (, for which one may close the contour in the lower half plane and 
take the contribution of the pole 0\ . 

This is reminiscent of the different regions encountered in the kinematics of the generalized parton 
distributions H(x, £, t), with the boundaries controlled by the relative values of x and £, where x and 
£ are analogous to our z and ( variables. 
The result of integration takes thus the form 



C J V DA Pi DA (h J Jo J DA Pi 

where DAp i denote the value of the product of propagators at the pole ft. 

A number of technical simplifications are now helpful. Firstly, one easily verifies that, at the leading 
logarithmic approximation we are interested in, the trace TrA\ may be simplified by taking the limit 
in 2 — > 0. Secondly, each integral over l 2 in Eq.® is UV divergent. However, it is a well-known 
classical result of QED that the sum of integrals corresponding to the six diagrams of Fig[T] is UV 
finite, so we separate UV divergent terms of each diagram in an algebraic way, and we show that the 
UV finiteness appears for the sum of diagrams A, B and C and for the sum of diagrams D, E and F 
separately. 

The traces are simple polynomials in l 2 , which may be written as Tr A\{0 — 0i) ~ aj 4 + jiQ 2 l 2 + 
SiQ 4 , where <Xi, and Si are dimensionless functions of z and £. Power counting in l 2 shows that the 
a^l 4 term in these integrals is ultraviolet divergent, since DAp i behaves as l 6 . Our aim is to recover 
the UV finiteness of the sum of diagrams before performing the integration over z. This cancellation 
of UV divergences occurs separately in each interval z £ [0, £] and z € [(, 1]. Let us concentrate for 
definiteness on the interval [0, £]■ The UV divergent part of I Al in this interval has the form : 

K r ji2 

-2iir dz J —af v (z,0, (8) 

where af m is a definite function of z and £. The contributions from diagrams B and C lead to 
similar equations as Eq.© with af m replaced respectively by bf m and cf lv . The cancellation of UV 
divergences in the sum of A, B and C diagrams is manifest through the fact that of v + bf v + cf v = 0. 
This relation allows us to extract the convergent part of the second term of Eq.® by adding and 
subtracting a term which reproduces the same UV divergence but which is IR finite: 

tM o- f „2 ( TrA 1 {0 = 1 ) l 4 dmf n \ 

I di l - dz J dl [—^ (1 2 +ro2) 3 °l M) 

- 2 ™l <dz J dl2 Q2T^ afV{Z > (9) 

The first term in Eq.Q is, by construction, both UV and IR finite and can be computed, whereas 
the second term in Eq.([9]) cancels out in the sum of diagrams A, B and C and thus does not need to 
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be evaluated. The same procedure is applied to the contribution coming from the interval z G [£, 1]. 
Moreover, the contributions from the diagrams D, E and F are considered in the same way. 

The final result of our calculation of the amplitude W\ reads, in the leading log Q 2 approximation 



Wi 



e 4 q Nc 

27T 2 ,„ 



dz [2z- 1) 



log- 



zC " • z( 

The amplitudes W2 and W3 are calculated in the same way and we get 



W 2 = 



(10) 



(11) 



and 



W 3 



e 4 q N C 
' 2tt 2 



dz 



^e{z-Q-U{c-z) 

zQ zQ 



-Mz-Q + ^6{l-z) 



log 



Q 2 



(12) 



In the following section we will interpret the results (fTU)) and (fT2"|) from the point of view of 
QCD factorization based on the operator product expansion, yet still in the zeroth order of the QCD 
coupling constant and in the leading logarithmic approximation. The crucial point is to note that 
the final contribution to this amplitude involves mixing of operators constructed from quark fields 
with operators constructed from photon fields [2]- This mixing can be understood by denoting the 
integrands of Eq. (fTU|) and Eq. (fT2"]) by J-(z, £) log ^p- and rewriting them by using the obvious identity 



Tfi 7TL ]\$ 

Hz,Q log^2 =F{z,Q logj^+T(z,C) log-^f 



(13) 



where Mp corresponds to an arbitrary QCD factorization scale. As will be shown below the first 

2 

term with log jjt may be identified with the quark GDA of the photon, whereas the second term 

with log -gf corresponds to the so-called photon GDA of the photon, coming from the matrix element 
of the correlator field from photonic fields which contributes at the same order in a em as the quark 
correlator to the scattering amplitude. The choice M F = Q 2 will allow to express the amplitude only 
in terms of the quark-antiquark fragmentation into two photons. 



3 QCD factorization of the DVCS amplitude on the photon 

We first consider two quark non local correlators on the light cone and their matrix elements between 
the vacuum and a diphoton state which define the diphoton GDAs $1, $3 : 

F q = J ^e i ^- 1 )i( 7 (p 1 )7(P2)|9(-|iV)7.iVg(|iV)|0) = \ 9 T^M)^M^x{z,Q^) (14) 

and 

F q = J ^e^-^hbxh^M-lAOT-A^ (15) 

where we note N — n/n.p, ^ v p n = e t±ua Pp a Np, with e 0123 = 1, and where we did not write explicitely, 
for simplicity of notation, neither the electromagnetic nor the gluonic Wilson lines. We will also define 
the matrix elements of photonic correlators 

F1 = j |H (2Z - I)f bfa)^)\F N *{-\N)F»{\N)\*) (16) 



4 



and 



pi 



dy 



27r . e ^-l)l( 7(pi)7(p2 )|^ ( _| 7V) ^ ( | Ar) | 0)) 



(17) 



where F Nft = N V F V ^ and F^ v = ±e^ p<T F pCT , which mix with correlators fH]) and (T5]). but contrarily 
to the quark correlator matrix element, are non zero at order oP em [2]. 

The quark correlator matrix elements, calculated in the lowest order of a em and as, suffer from 
ultraviolet divergences, which we regulate through the usual dimensional regularization procedure, 
with d = 4 + 2e. We obtain (with I — - + 7B - log47r) 



Nce 2 q 
4tt 2 



9t" e *u(Pl) e t(P2) 



— + log TO Z 

e 



(18) 



with 



= - C) + ^^^(^ - C) + ^^^^(C - + ^^^^(C - (19) 



c 



c 



for the /_t <-» v symmetric (polarization averaged) part. The corresponding results for the antisymmetric 
(polarized) part read 



F" 



Air 2 



log 77?/ 



with 



F(z, = j6{z - - j6(z ~ C) ~ jO(C -z) + jO(C ~ z). 



(20) 



(21) 



Let us stress again that here we concentrate only on the leading logarithmic behaviour and thus focus 
on the divergent parts and their associated logarithmic functions. The ultraviolet divergent parts 
are removed through the renormalization procedure (see for example [5]) involving quark and photon 
correlators (O q , O 7 ) corresponding to one of the two following pairs 



{q{- V -N) 1 .Nq{lN), F N »{- V ~N)F% {\n)) , 
(g(-|7V) 7 .7V 7 5 g (|iV), F N ^- V -N)F^{ V -N)). 



(22) 
(23) 



The renormalized operators are defined as: 



O" 

Of 



O" 
0~ 



The matrix element of the renormalized quark-quark correlator is thus equal to 

< 7(Pi)7(P2)|O||0 >= Z qq < 7 (pib(P2)|O g |0 > +Z qi < 7 ( Pl ) 7 (p 2 )|O 7 |0 > 



with Z qq = 1 



O 



(*)• 



(24) 



(25) 



Since the matrix clement < 7(^1)7(^2) |O 9 |0 > contains a UV divergence 

(see Eas.lfTB)). (|2U)) ) and since < 7(^1)7(^2) |O 7 |0 > is UV finite and of order a° m , one can absorb this 
divergence into the renormalization constant Z qi . The normalization of the renormalized correlator 
is fixed with the help of the renormalization condition which is chosen as 



<7(Pi)7(P2)|O|j|0>=0 at M R = m. 
In this way the renormalized GDA with vector correlator is equal to 



n 



Ncel 



in 



-^L g ^e;( Pl K(p 2 ) log j^F[z, 



(26) 



(27) 
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and wc have a similar result for the renormalizcd GDA with axial correlator 



F R = -^{-i^YMYMlog^FizX) ■ (28) 

As we want to use QCD factorization formula, which correspond to the factorization scale Mp, we 
identify now the renormalization scale Mr with Mp, 

M R = M F . (29) 

Eqs. (|2"Tll2"5|) together with Eqs (fT§ll2"Tj) permit us to write the expressions of the diphoton generalized 
distribution amplitudes: 

Nn € 2 m 2 

$?(z,C,0) = —g-JL]o g ^F(z,Q, (30) 

F 

$l(z,C,0) = -- ^llog^F^C). (31) 

F 

Now we are able to write the quark contribution to the 7*7 — > 77 amplitude at threshold as a 
convolution of coefficient functions and GDAs &f(z, £, 0) 

1 1 



W* = / dzC«,(z)$l(z,(,0) , W§ = dzC q A (z)<t> q 3 (z,(,0) , (32) 





where the Born order coefficient functions Cy, A attached to the quark-antiquark symmetric and 
antisymmetric correlators are equal to: 

We recover in that way the In j^- term in the right hand side of Eq. (fTlff . The photon operator con- 
tribution involves a new coefficient function of order calculated at the factorization scale Mp, 
which plays the role of the infrared cutoff. This function is convoluted with the photonic correlators 
F N n(-%N)F*(%N) and F N ^(-^N)F^ (fiV) ( c f Eqs.fHi [T7J)) which are of order a° em . These con- 
volutions have the same expressions as those in Section 2 (Eqs.lfTUl fT2"])) with the quark mass replaced 
by the factorization scale, m — * Mp. In this way, we recover the second term in the right hand side 
of Eq. (fl"3|) . This equation in fact reflects the independence of the scattering amplitude on the choice 
of the scale Mp, which is controlled by the renormalization group equation. 

The factorization scale Mp can be chosen in any convenient way. The choice Mp = Q 2 kills the 
logarithmic terms coming from the contribution of photonic GDAs, so that the scattering amplitude 
is written (at least in the leading logarithmic approximation) solely in terms of the quark correlators. 
With Mp = Q 2 we interpret this process within the parton model (see e.g. [10], [3] for an analogous 
interpretation of parton distributions inside the photon). 

4 The diphoton GDAs and their QCD evolution equations 

We have thus demonstrated that it is legitimate to define the Born order diphoton GDAs at zero W 
and at Mp — Q as 

*; ( , >f ,o) = ^io,^ 2 ^-0^ ^,^,-0 



2tt 2 b m 2 



g(*-0+ , J 0{z-Q 

(34) 



G 



and 



C,o) 



Nee 2 g 2 
2ir z 



^9(C-*) + j6(C-z) 



(35) 



Since we focus on the leading logarithmic contribution, we only obtain the anomalous part of these 
GDAs. Their z— and (— dependence are shown on FigsJDandEB 




Figure 2: The unpolarized anomalous diphoton GDA <I>^ 2-k 2 /{Nc e 2 log(Q 2 /m 2 )) at Born order and at 
threshold for ( = 0.1 (dashed), 0.2 (dash-dotted), 0.4 (solid). 




Figure 3: The polarized anomalous diphoton GDA $3 2-k 2 /{Nc e 2 log(Q 2 /m 2 )) at Born order and for ( = 0.1 
(dashed), 0.2 (dash-dotted), 0.4 (solid). 



Note that Qj{z,(,0) GDAs are discontinuous functions of z at the points z = £ and z = £. 
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Nevertheless, these GDAs still verify the property of polynomiality 

n+l 



dz (2* -1)™ <!>?(*, C) = J2 a ^ k > 

„1 n+l 

/ dz (2z-irm z , o = E a ^ fe - 

Ja k=0 



(36) 
(37) 



As we have two photons in the final state, the only non vanishing correlations are the C-even i.e. singlet 
sector of the GDAs which correspond to the combinations of operators ^(O q (xi, X2) — O q (x2, xi)) 
and i(0 9 (xi, x^) + O q (x2, xi)). This singlet sector means for the GDAs 5>f the combinations 

$«(*, C, 0) = ±($?(z, C, 0) - C, 0)) and C, 0) = |(3|(z, C, 0) + , C, 0)). From 

the point of view of flavor decomposition, we can distinguish either singlet sector 



Nf N f 
9=1 <J=1 



or non-singlet sector 



(38) 



(39) 



where q and g' are two different quark flavors. 

Let us illustrate the effect of evolution in the simplest, i.e. without mixing with gluons, case of 
non-singlet and vector GDA <f>^ s . 

Switching on QCD, the non-singlet and vector sector of diphoton GDAs evolves according to the 
ERBL evolution equation [H] modified by the presence of the anomalous part 



Q 2 ^<S>'? S (zX,Q 2 )= £ duV NS (u,z,Q 2 )<S>>? s (u,(,Q 2 ) + (e 2 q ~e 2 q ,)f[(z,(), (40) 

with $^ s (z,C,<2 2 ) = z(l-z)$f s (z,C,Q 2 ) and/i(z,C) = z(l- z)f[(z,(). Here /1 C) is defined by 
the r.h.s. of Eqs. ([3~4"|) , as the corresponding functions which multiply e 2 In Q 2 /m 2 . The QCD kernel 
Vns is known [3] and at the leading order accuracy it reads 



V NS (u,z,Q ) = — C F 

Z7T 



-9{z - u) 1 + 



-9{u-z) 1 



(41) 



with Cf = (Nq — l)/(2Nc) and the + prescription is 



da 



/(<*) = 



da 



f(a) - /(0) 



(42) 



The kernel Vats is diagonalized in the basis spanned by Gegenbauer polynomials Cp 3 ^ (x) [TT] . This 
property permits us to write the solution of the equation (|4TJ|) as 



00 

p odd 



(43) 



A 1, 



log — o 



c -lqq(P) 



- ,2 2 ,, q 2 m 

e Q - e a , ) log — ^ rT - 

1 ' u 33-2Af f 



C( 3 / 2 '(2z-l), 



where ^4 p are integration constants, TV/ is the number of flavors, 7 qg (p) are the usual anomalous 
dimensions 

(\ 1 p+1 



(44) 



fc=2 
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and the coefficients f p (() are projection of f[(z, £) on appropriate Cp 3 ^ 2 ^ polynomials 

W = {p+ 2 l)t + 2) l!^ ^^OC^^-D- (45) 

Since the expressions for f' p (C) are lenghty, we do not presente them explicitely. For large Q 2 , the 
solution (|43|) . at the leading logarithm level, reads 

C,Q 2 )^ (eg -4) log g zz £ 1 /^ (p) C^(2z-1). (46) 

p odd 1 + |3 33-2JV / 

The formula (|46|) shows the known result that the anomalous part of GDAs dominates at large Q 2 . 
The z, £ dependences of $^ are nevertheless modified by strong interaction, due to the presence of the 
denominator 1 + 6 • 

The result of this method, for Nf — 2, is shown in FigfJJ The procedure of decomposing a GDA 
into series of Gegenbauer polynomials ensures a good numerical description of both the unevolved 
GDA and the one with the evolution taken into account, except in the regions z close to £ and 1 — £ 
(it is well known that theta function are very badly described by series of polynomials) . Fig|4] was 
obtained after combining the sum of 41 (201) contributions of Gegenbauer polynomials for the regions 
< z < £ and 1 — £<z<l(£<z<l — £). Such a choice for the regions < z < £ and 1 — £ < z < 1 
is motivated by numerical instabilities of the Gegenbauer series (|4"B"|) when z is close to and 1, which 
requires to truncate the series at a moderate number of terms (41 terms), which nevertheless leads 
to a good stability of results. These numerical instabilities are due to the fact that, as soon as z 
deviates from or 1, f p (Q (Eq. HS]) ) becomes very large with £ fixed and p larger than 50. This is 
in contrast with the region £ < z < 1 — ( where a very good stability of the results is achieved when 
summing a very high number of terms (we took 201 terms: taking more terms would then generate 
huge instablities in the vicinity of £ and 1 — £) . 

As can be seen from FigJH QCD evolution affects the form of the GDA in the whole z range. 




Figure 4: The bare (solid line) and the QCD evolved (dashed line) non-singlet and vector sector of diphoton 
GDA at large Q 2 $^ s /((e^ - e\) \og{Q 2 /m 2 )) (see Eq.®}) for £ = 0.2. The dotted line is the truncated (see 
text) Gegenbauer expansion of the bare GDA. 



In practice, one cannot trust the truncated Gegenbauer expansion near the discontinuity points. 
To obtain a trustable behaviour of the evolved GDA around z = £ and z = £, one has to sum in 
Eq. (|46p the infinite series involving Gegenbauer polynomials, which is a non trivial task. Because of 
that, one could think about a method based on direct iteration of the non-homogeneous equation (|40p . 
keeping only the leading terms at large Q 2 and we refer for all details to the Appendix. In particular, 
we have solved the evolution equation (|40l I41| in the vicinity of singularity points z = £, 1 — £. Because 
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of the antisymmetry of the GDA let us concentrate our discussion on the form of the evolved 

GDA around the singular point z = £. The iteration of the kernel (|40[) generates leading terms of the 
form K n ln™(C — z) for z — > £~ and K n ln™(z — £) for z — » £ + > which we resummecO] in 



*r(*-r,0 = 5 (e -<)log^ 



4C 2 - C - i 



(l-Xln(C-z)) + 



and 



1 ^ 



4C 2 - C - i 

c 



(l-Xln(z-C)) 



l-ifln(C-z) 



1 



(47) 



(48) 



where K = 6 Cp/(UNc—2Nf) . We display in FigEJthe asymptotic resummation given by Eqs. (|47|48[) . 
As mentionned above, the resummation of logarithmic contribution is required in the regions K log \z — 



0.399 0.3995 



-0 . 5 



-1 



-1 . 5 

-2 
-2/5 



.4005 0^401 



Figure 5: The behaviour of the evolved GDA from the leading logarithmic resummations (|47l 148 [I in the 
vicinity of the discontinuity point z = £, for the case f = .4 (dotted curve). The Born GDA is given for 
reference as the solid line. 

C| > 1 and if log \z-(\> 1, which in our case (K = 8/29) corresponds to \z - £\,\z - C| < 3._1(T 2 . 
In practice, subleading contributions are non negligible in a larger domain in z around C and and 
prevent us from getting trustable numerical results based on this iterative method. This is due to 
the numerically large values of K. This iterative method stabilizes only when applied for unphysically 
small values of K (typically < .03). 

Note that these singularities of the GDA at z = ( and z — 1 — C does not affect the DVCS 
amplitude (f3"2"|) since they are integrable singularities!!. 



Conclusions 



We derived the leading amplitude of the DVCS (polarization averaged or polarized) process on a 
photon target at threshold. We have shown that the amplitude coefficients Wf factorize in the forms 
shown in Eq. (|32p . irrespectively of the fact that the handbag diagram interpretation appears only 
after cancellation of UV divergencies in the scattering amplitude. We have shown that the objects 
&1{z, CjO) are matrix elements of nondocal quark operators on the light cone, and that they have 
an anomalous component which is proportional to log(Q 2 /m 2 ). They thus have all the properties 
attached to generalized distribution amplitudes, and they obey non-homogeneous ERBL evolution 
equations. This new type of evolution equations is an interesting playground to study the effects of 
gluon radiation on a non diagonal object such as a GDA. 



^^Note that this series corresponds to the resummation of soft gluons exchanged between quark-antiquark lines, when 
iterating the ERBL kernel on the bare GDA. 

2 The coefficient function Cy l|33| l is regular. This would remain true after taking into account its QCD evolution 
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Appendix 

Using 

'■ (9i) ^.»c-"^r 1491 

denoting t = log and redefining the GDA as $ q (z, £, t) — zz & q (z, £, t), the evolution equation (|4"0"|) 
reads 

^< s (z, C, t) = y cfc ^(u, z) < 5 (", C, t) + (e^ - 4)#(*, , (50) 
where f[(z, Q, with the help of Eq. ([34"]) . is 



ft r\ Nc 

Z7T Z 



(51) 



The solution of Eq. (j50|) is the sum 

&» s (z,(,t)=g(z,<;,t) + i(z,(,t), (52) 
where g(z, t) is the general solution of the homogeneous evolution equation, 

d f 1 

—g(zX,t) = J duV NS {u,z,t)g{u,C,,t) (53) 
and i(z, £, t) = t h{z, Q is a particular solution of the non-homogeneous evolution equation 

h(z,0=t [ duV NS (u,z)h(uX) + h {0 \z,0, (54) 



where h^(z,() — (e 2 , — e \')fi( z ,Q • As we w iU De interested in the leading terms for large Q 2 , the 
general solution of the homogeneous equation (|53")) can be omitted. Eg. (f5"4")) reads 

h(z,0 = K [ duV(u,z)h(u,() + h (0) (z,(), (55) 







with if = 6 Cf/(11Nc — 2Nf) and V is the expression inside [■ ■ • ] in Eg . (|4"T]) . This equation may be 
solved numerically by the method of successive iterations, which represents the solution as 

oo 

h( Z ,Q = J2( k vt ® h(0) ■ ( 56 ) 

n=0 

Since we discuss the vector GDA, satisfying Eq. (|50")) . which is antisymmetric with respect to the 
substitution z — ► z, it is sufficient to construct the solution for z € [0, 1/2] . Similarly, the solution is 
symmetric with respect to the substitution £ — > £, and thus without loss of generality, we restrict our 
analysis to £ €E [0, 1/2] . 

We write the solution as 

h(z, = h^z, 0(C - 2) + ft 2 (*> 0(z - C) , (57) 
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and use the same decomposition for h^°'(z, () 

(z, C) = (z, C) 0(C -z) + 4 0) (z, C) 9{z - C) , 
where, with the help of Eq. ([5T]) . 

-2] and 4 0) (z,C) = (^-4) 



^c) = (^-4)i(^ 

In these notations, Eg. ([55]) reads 



C (gu - 1) 

£ MM 



/H(z,0 



A" 



1 - 2z Ml 
it 



z z — u 
l-2u u 1 



f/»i(«,0-fti(^C) 



z 



+ 



z 
A" 



1/2 r 



1 - 2u 



z z — u 
u 



z — u 



du 



(58) 



(59) 



(60) 



du 



u — z z — u 



1 -2z 



z z 
1 - 2z 



-u + u 



h 2 (u, Qdu 
1 



hi(u, Q)du 



zz 



uh 2 {u,Q) + 



z(z — u) z{z — u) 
*h 2 (u,0-h 2 (z,0 *§h 2 (u,0+h2(z,0 



(61) 



1/2 r 



1 -2u 



h 2 (u,() 



z — u z — u 

?§h 2 (u,() + h 2 (z,0 , ffcaKO-fo^O 



In order to understand the nature of the singularity observed at z = £ and z = it is useful to 
observe that the dominant contribution appearing at each iteration involves K n ln™(£ — z) for z — > £~ 
and A n ln"(z — £) for z — > £ + . Our aim is now to sum all such leading logarithmic terms. For that, 
we note that such terms are generated by the approximate evolution equations, obtained from (|60p 
and (|6ip when keeping only terms generating large In™ \z — C| '■ 



ht(z,0 = K 
~h 2 (z,0 = K 
where h£ = (2( - 1)/C 2 and h 



hi(u,Q - hi(z,Q 



du 



1/2 



h 2 (u,Q 



du 



C fei(«,C) , f 1/2 h 2 (u,0-h2(z,0 



+ 



z — u 



du 



+ h 
+ h+, 



q — Hq — 1/ (CC) 1 see Eq. ([5^1 . The solution of these equations are 



&i(*>0 - ~ ((V + fe + )(i - KH( - z)) + i _^ ln( ^°_ ^ 



1 



h Q -h+ 



1 - Aln(z-C) 



Identifying 



$^ S (z->r,C) = zz^-e^log^ft^z.O. 



^ s (z^C+,() = zz(e 2 q -e 2 q ,)log^h 2 (z,C) 



(62) 
(63) 

(64) 
(65) 

(66) 
(67) 



we recover Eqs. (|4"?I |4"5|) of section [U 

Before proving that (|6"4"1 (HSJ) are required solutions, let us explain first, in a somehow heuristic 
way, how one can guess the form of solutions ([64"! l65|) from Eq. ([62|) and (|63|) . We rewrite the second 
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(resp. first) integral of Eq.([H2]) (resp. (JB3J)) using h 2 (i) ( u , C) = ^2(1) ( u > - ^2(1) ( z , C) + ^2(1) (z> C) and 
perform the integration of the constant term at leading logarithmic order, which gives 



ht(z,C) =K 
h 2 (z,() =K 



-du 



thxM-hi&Qj , /" 1/2 h 3 (u,C)- jfcfoC) 



-dit 



u — z 



Kh 2 (z,C) ln(C - z) + V , (68) 
Kh^ZiQ HC-z) + h+. (69) 



The trick, valid at leading logarithmic order, is now to replace in the denominator of the integrand 
of the right hand side of these equations z by £ and to replace the lower (resp. upper) bound of the 
second (resp. first) integral of Eq. lpJS]) (resp. (f6"9"|) ) by z, which gives 



h 1 (z,() = K 
h 2 (z,()=K 



/Li(«,C)-/ii(z,C) du , f 1/2 ^,0-M^C) 



u-C 



-du 



hi(u,Q ~ hi(z,() t f 1/2 h 2 (u,0 ~ h 2 {z,() 



C-u 



u-C 



Writing symbolically 



£ = /ii + h 2 and A = hi — h 2 , 
these two equations decouple in leading logarithmic accuracy into 

E(«,C)-E(«,C) 



-tf&afoO ln(C - 2) + />o , (70) 
-tffei(«,C) ln(z - C) + h+ . (71) 

(72) 



E(z,C) = 2K 
The algebraic equation (f?3"|) is solved in 

A(z,C) = 



C-u 



du- K ]n(z -() + h Q + 



1 - K ln(z - C) 



(73) 
(74) 

(75) 



Eq. (fT4f can be turned into an elementary differential equation when integrating the second term of 
the integral and differentiating with respect to z, giving 



Yf\l-K ln(C - z)] + E(z) = 2K 



which solution is 
From Eq. (fT2"j) . 



C-z — C-z' 
E(z,C) = C[l-Xln(C-z)]. 



hi = — i — and h 2 = — - — 



(76) 

(77) 
(78) 



and the constant C is fixed as C = h + by comparison of the lowest order term in the K expansion 
with Eqs.dHlEU). This justifies the ansatz ([Ml IBS]) . 

We now prove that the expressions (|6"4"1 1^5]) are the wanted solutions, by direct substitution of 
them into the equations (|62| I63p and performing necessary integrations. Let us consider as example 
the Eq.([nil) for h%(z, Q. The right hand side of it takes thus the form 



K 2 

rhs of Eq. dS2J) = — 



du 



-(K + h+) In 



(h Q -h+) InfEf 



K 

y 



1/2 



du 



u — z 



[ho +h+)(l -K ln(u-C))- j 



C-z [1 - K ln(C-u)][l- K ln(C-z) 



JCln(u-C) 



(79) 
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The three integrals which appear in (|79p are calculated in a straightforward way 

-Li 2 (l 



du C — u 
In ■ 



z — u ( — z 



s " :) — ^ 2 (C-*), 



C — z c-"-*o 2 



(80) 



1/2 



, 1/2 -z 
= In ■ 



( — z c-«^o 



-ln(C-z), 



(81) 



1/2 



— In(u-C)=In(V2-0Ml + ^— ~) + Li2 / ' 1/2 C 



£ — z J c-«^o 



-l/21n 2 (C-z). (82) 



1/2 



The integral J i-K\n{u-Q can ^ e eva l ua ted by developing i_k\^{u-q * n ^ s Taylor series in 

K\n(u — C) and by performing integration over u in each term of this expansion. The leading terms 
in the limit of small values of £ — z are then summed up again, which leads to the formula 



1/2 



du 



1 



1 



- z 1 - K\n{u - C) C-«^o K 



In [1 - K ln(C - z) 



(83) 



More involved is the evaluation of the integral 

z 

, , du 



C-3 



z — u 1 — K ln(f — u) 



We calculate it by performing first the change of variables t = y^- which leads to the expression 



K 



dr 



In - 



IMC-*) 



r(l -r) 1 + Glnr 



with G 



K 



l-K\n{Q-z) 



(84) 



Next, we represent the integrand of / as a double series in r and Glnr and perform the integration 
over r in each term of this double series by keeping leading terms in the limit Q — z — > : 



G 
K 



n=0 p=0, . n=0 



P=0 ^(c- Z ) 



ln" +J (C-z) 
n + 2 



K 2 



(l-A-ln(C-z)) 



ln(l -ifln(C - z)) + 



ifln(C - z) 



1- A'ln(C-z) 



(85) 



By substituting expressions ((501 IST1 1521 [551 155)) into Eq. (j?9"|) we recover the expression (jo"4")) for /ii(z, £). 
In a similar way we have checked that the expressions (|rH|) and (fBl)]) satisfy the second approximated 
evolution equation (f63|) for /i2(z,£). 
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